PATENT APPLICATION 
REAL-TIME MEASUREMENT OF ULTRASHORT LASER PULSES 

BACKGROUND OF THE INVENTION 

1. Field of the invention 

[0001] The present invention relates to a method and apparatus for the 

real-time measurement of ultrashort laser pulses based on frequency- 
resolved optical gating (FROG). 

2. Background Art 

[0002] Ultrashort optical pulses have very short time durations, 

typically less than a few 10's of picoseconds. As a result, ultrashort optical 
pulses are spectrally broad. Because the index of refraction of materials is a 
function of wavelength, different wavelengths of light travel at different 
speeds in optical materials causing the properties of ultrashort optical pulses 
to change as they propagate. The shape of the pulse can influence how the 
pulse itself interacts with materials, further complicating the analysis of 
ultrashort optical pulse transmission and propagation. Ultrafast laser pulse 
measurement tries to obtain both the intensity profile of the pulse as well as 
the "phase of the pulse," which is the actual variation of the frequencies that 
make up the pulse. 

[0003] Generally, time resolved measurements of events use shorter 

events to resolve variations in the characteristics of the events as a function 
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of time. Unfortunately for the analysis of ultrashort optical pulses, shorter 
events do not exist and modern electronics are insufficiently fast to allow for 
direct time resolved measurements. One effective measurement technique 
for ultrashort optical pulses is called "frequency-resolved optical gating 7 ' or 
FROG. FROG is based on the observation that an optically formed 
spectrogram contains all of the information about the pulse to be measured. 
A two-dimensional phase retrieval process extracts the pulse from its 
spectrogram. The main disadvantage of frequency-resolved optical gating is 
that the pulse extraction using phase retrieval can be slow, which is 
undesirable for real-time applications. Aspects of the present invention allow 
FROG to be implemented as a real-time pulse measurement technique that is 
both accurate and fast. 

[0004] Mathematical Representation of an Optical Pulse 

[0005] The mathematical representation of ultrashort optical pulses is 

discussed to provide background for the discussion of the FROG technique. 
The time-dependent variations of an optical pulse are embodied in the pulse's 
electric field, A(t), which can be written' 

A(t) = Re[E(t)e i6,ot ] (1) 

where coo is the carrier frequency and Re refers to the real part. A(t) can be 
used in this form for calculations, but it is generally easier to work with a 

different representation that removes the rapidly varying coo part, e'^ 0 ' , and 
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uses as a representation a slowly varying envelope together with a phase 
term that contains only the frequency variations. This representation, which 
does not include the rapidly varying carrier frequency, is 

E(t) = [I(t)] m e i<p{0 (2) 

where Kt) and <p(t) are the time-dependent intensity and phase of the pulse. 
(EXx) is complex.) The frequency variation, Q(t), is the derivative of cp(t) with 
respect to time- 

Q(t) = -d<p(t)/dt (3) 

[0006] The pulse field can be written equally well in the frequency 

domain by taking the Fourier transform of equation 2- 

E(o)) = [I(a)] U2 e-'* {a,) (4) 

where Ku>) is the spectrum of the pulse and cpfo) is its phase in the frequency 
domain. The spectral phase contains time versus frequency information. 
That is, the derivative of the spectral phase with respect to frequency yields 
the time arrival of the frequency, or the group delay. 

[0007] Obtaining the intensity and phase, Xt) and <p(t) (or j(co) and cp(co)) 

is called full characterization of the pulse. Common phase distortions include 
linear chirping, where the phase (either in the time domain or frequency 
domain) is parabolic. When the frequency increases with time, the pulse is 
said to have positive linear chirp; negative linear chirp is when the high 
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frequencies lead the lower frequencies. Higher order chirps are common, but 
for these, differentiation between spectral and temporal chirp is required 
because spectral phase and temporal phase are not interchangeable. 
[0008] Frequency Resolved Optical Gating 

[0009] Frequency-resolved optical gating (FROG) measures the 

spectrum of a particular temporal component of an optical pulse, shown in 
FIG. 1, by spectrally resolving the signal pulse in an autocorrelation- type 
experiment using an instantaneously responding nonlinear medium. As 
shown in FIG. 1, FROG involves splitting a pulse and then overlapping the 
two resulting pulses in an instantaneously responding x (3) or x (2) medium. 
Any medium that provides an instantaneous nonlinear interaction may be 
used to implement FROG. Perhaps the most intuitive is a medium and 
configuration that provides polarization gating. 
[0010] For a typical polarization- gating configuration, induced 

birefringence due to the electronic Kerr effect is the nonlinear-optical 
interaction. The "gate" pulse causes the x (3) medium, which is placed 
between two crossed polarizers, to become slightly birefringent. The 
polarization of the "gated" probe pulse (which is cleaned up by passing 
through the first polarizer) is rotated slightly by the induced birefringence 
allowing some of the "gated" pulse to leak through the second polarizer. This 
is referred to as the signal Because most of the signal emanates from the 
region of temporal overlap between the gate pulse and the probe pulse, the 
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signal pulse contains the frequencies of the "gated" probe pulse within this 
overlap region. The signal is then spectrally resolved, and the signal 
intensity is measured as a function of wavelength and delay time x. The 
resulting trace of intensity versus delay and frequency is a spectrogram, a 
time- and frequency-resolved transform that intuitively displays the time- 
dependent spectral information of a waveform. 
[0011] The resulting spectrogram can be expressed as- 



where E(t) is the measured pulse's electric field, g(t-x) is the variable -delay 
gate pulse, and the subscript E on Se indicates the spectrogram's dependence 
on E(t). The gate pulse g(t) is usually somewhat shorter in length than the 
pulse to be measured, but not infinitely short. This is an important point- an 
infinitely short gate pulse yields only the intensity I(t) and conversely a CW 
gate yields only the spectrum I(co). On the other hand, a finite -length gate 
pulse yields the spectrum of all of the finite pulse segments with duration 
equal to that of the gate. While the phase information remains lacking in 
each of these short-time spectra, having spectra of an infinitely large set of 
pulse segments compensates this deficiency. The spectrogram nearly 
uniquely determines both the intensity I(t) and phase cp(t) of the pulse, even if 
the gate pulse is longer than the pulse to be measured. If the gate is too long, 
sensitivity to noise and other practical problems arise. 




(9) 
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[0012] In FROG, when using optically induced birefringence as the 

nonlinear effect, the signal pulse is given by" 



E sig (t,T)ccE(t)\E(t-T)\ 2 ( 10 ) 

So the measured signal intensity Ifrog(©,x), after the spectrometer is" 

r) = | £ E(t)\E(t - rfe- im di\ (1 1) 

The FROG trace is thus a spectrogram of the pulse Ekt) although the gate, 
| Efo'x) | 2 , is a function of the pulse itself. 

[0013] FROG is not limited to the optical Kerr effect. Second harmonic 

generation (SHG) FROG can be constructed to analyze relatively weak pulses 
from oscillators and is typically more sensitive than polarization- gating 
FROG. For SHG FROG, the pulse is combined with a replica of itself in an 
SHG crystal as illustrated in FIG. 2. 

[0014] To see that the FROG trace essentially uniquely determines Efct) 

for an arbitrary pulse, note that Mt) is easily obtained from Esi/t, x). 
Equation (ll) can then be written in terms of E S iJ<t, Q), the Fourier transform 
of the signal field E S i^t y x) with respect to delay variable x. This gives the 
following, apparently more complex, expression* 
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Equation (12) indicates that the problem of inverting the FROG trace 
Ifrog((o,i) to find the desired quantity Esi^t, x) is that of inverting the 
squared magnitude of the two-dimensional (2-D) Fourier transform of 
Esi^ty x). This problem, which is called the 2-D phase -retrieval problem, is 
well known in many fields, especially in astronomy, where the squared 
magnitude of the Fourier transform of a 2-D image is often measured. At 
first glance, this problem appears unsolvable; after all, much information is 
lost when the magnitude is taken. The 1-D phase retrieval problem is known 
to be unsolvable (for example, infinitely many pulse fields give rise to the 
same spectrum). Intuition fails badly in this case, however. Two- and 
higher-dimension phase retrieval processes essentially always yield unique 
results. 

[0015] FROG Inversion 

[0016] An iterative 2-D phase retrieval process is used to extract the 
pulse information from the measured FROG trace as illustrated generally in 
FIG. 3. This phase retrieval process converges to a pulse that minimizes the 
difference between the measured and the calculated FROG trace. 
Application of this phase retrieval process to FROG has been problematic in 
the past. Some recent applications use a generalized projections algorithm, 
which converges quickly, along with faster computers, to track pulse changes 
at rates of 20 Hz or greater, making FROG a real-time pulse measurement 
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technique. Indeed, programs for analyzing FROG traces are commercially 
available. 

[0017] The original FROG inversion process, using what is commonly 

referred to as the vanilla algorithm, is simple and iterates quickly. On the 
other hand, the process tends to stagnate and give erroneous results, 
especially for geometries that use a complex gate function such as SHG or 
self-diffraction. Improved strategies using different algorithms, including 
brute force minimization, were developed to avoid stagnation, at the expense 
of both iteration speed and convergence time. Later a numerical technique 
called generalized projections brought a significant advance in both speed 
and stability. The generalized projections technique proceeds after an 
iteration by constructing a projection that minimizes the error (distance) 
between the FROG electric field, Esijti, x), obtained immediately after the 
application of the intensity constraint, and the FROG electric field calculated 
from the mathematical form constraint. The projection constructed in this 
manner is used as the starting point of the next iteration. 
[0018] The first implementations of the generalized projections 

technique used a standard minimization procedure to find the electric field 
for the next iteration (which can still be slow). For the most common FROG 
geometries, PG and SHG, there are substantial advantages to a different 
strategy that directly determines the starting point for the next iteration. 
This strategy, called Principal Components Generalized Projections (PCGP), 



8 



converts the generalized projections technique into an eigenvector problem. 
The PCGP technique has achieved pulse characterization rates of 20 Hz. 
[0019] The goal of phase retrieval is to find the Mt) that satisfies two 

constraints. The first constraint is the FROG trace itself which is the 



magnitude squared of the ID Fourier transform of E S i g (t,x)- 



FROG 



.o-|r. 



E sig (t,T)e- ia "dt 



(13) 



The other constraint is the mathematical form of the signal field, E S i£x,x), for 
the nonlinear interaction used. The mathematical forms for a variety of 
FROG beam geometries are: 



E sig (t,T)cc 



E(t)\E(t-rf 

E(t) 2 E\t-r) 

E(t)E(t-r) 

E(t) 2 E(t-r) 



PG FROG 

SD FROG 
SHG FROG 
THG FROG 



(14) 



where PG is polarization gate, SD is self-diffraction, SHG is second harmonic 
generation and THG is third harmonic generation FROG. 
[0020] All FROG implementations work by iterating between two 
different data sets: the set of all signal fields that satisfy the data constraint, 
iFRodss, x), and the set of all signal fields that satisfy equation 14. The 
difference between FROG implementations is how the iteration between the 
two sets is completed. In the case of generalized projections, the .E(t)'s are 
chosen such that the distance between the £kt) on the magnitude set and the 
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Mt) on the mathematical form set is minimized. This is accomplished by 
minimizing the following equation: 

where E^ {DC) (t l9 tj) is the signal field generated by the data constraint, 

and£^ + { ^ F) (r / ? r y ) is the signal field produced from one of the beam geometry 

equations 14. For the normal generalized projections technique, the 
minimization proceeds using a standard steepest decent algorithm; the 
derivative of Zwith respect to the signal field is computed to determine the 
direction of the minimum. The computation of the derivatives is tedious! the 
derivatives are tabulated in, for example, Trebino, etal, Rev. Sci Instrum., 
68, p. 3277 (1997). 

[0021] An alternative to the minimization of equation 15 is principal 

components generalized projections (PCGP). PCGP computes the starting 
point of the next iteration through an eigenvector problem, reducing the 
computation for the next iteration to simple matrix-vector multiplies. PCGP 
works for both the PG and SHG beam geometries, is simple to program and is 
fast, as described in D. J. Kane, IEEE J. Quant Elec, (1999). 
[0022] Self Checks in FROG Measurements 

[0023] Unlike other pulse measurement techniques, FROG can provide 

a great deal of feedback about both the quality of the measurement 
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(systematic errors) and the quality of the technique's performance. A good 
check for convergence is the FROG trace error together with a visual 
comparison between the retrieved FROG trace and the measured FROG 
trace. The FROG trace error is given by- 



where a is a renormalization constant, Ifrog is the measured FROG trace 
and Irfiocftico, t) is computed from the retrieved electric field. Typically, the 
FROG trace error of a PG measurement should be less than 2% for a 64 x 64 
pixel trace, while the FROG trace error of a 64 x 64 pixel SHG FROG trace 
should be about 1% or less. Acceptable FROG trace errors decrease as FROG 
trace size increases and increase for smaller FROG traces. These values are 
only guidelines and variations are expected in specific circumstances. For 
example, acceptable phase retrievals of large and very complicated FROG 
traces can produce larger FROG trace errors. 

[0024] In a good FROG measurement, the spectrum of the retrieved 

pulse should faithfully reproduce the salient features of the pulse's measured 
spectrum. SHG FROG even provides an additional check called the 
frequency marginal. The sum of an SHG FROG trace along the time axis 
yields the autoconvolution of the pulse's spectrum providing two ways the 
FROG measurement can be checked. First, the autoconvolution of an 




(16) 
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independently measured spectrum can be compared to the sum of the FROG 
trace along the time axis (the frequency marginal) providing an indication of 
how well the measurement was made. For example, if the doubling crystal 
was too thick in the pulse measurement, the FROG trace's frequency 
marginal will be narrower than the autoconvolution of the measured 
spectrum. Second, comparing the autoconvolution of the retrieval pulse 
spectrum with the FROG trace marginal can provide a test of convergence in 
addition to a test of the measurement. Phase matching problems show up as 
a mismatch between the FROG trace marginal and the autoconvolution of the 
retrieved spectrum. 

[0025] Because FROG is a spectrally resolved autocorrelation, 

summing any FROG trace along the frequency axis yields the autocorrelation 
of the measured pulse. This autocorrelation can be compared to an 
independently measured autocorrelation, or a comparison can be made 
between the frequency sum of the FROG trace and the autocorrelation 
calculated from the retrieved pulse to determine algorithm convergence and 
the quality of the measurement. 
[0026] Making FROG Real-Time 

[0027] In order to make FROG a real-time pulse measurement 

technique, the FROG device, the data acquisition, the 2-D phase retrieval 
process, and the user interface are integrated. One real-time ultrashort laser 
pulse measurement implementation combined a multishot FROG device with 



12 



two digital signal processing (DSP) boards in a fast personal computer. This 
implementation is described in D. J. Kane, IEEE J. Quant EJec, (1999). One 
DSP card was devoted to data acquisition and the other DSP card was used 
to make the calculations. 

[0028] The DSP implementation was considered expensive for some 

applications and so efforts were directed to develop a FROG-based ultrashort 
laser pulse measurement that did not use DSP cards. One implementation 
that did not use a DSP, called VideoFROG, used a standard CCD video 
camera to acquire the FROG trace and a standard computer frame grabber to 
digitize the video signal. The computer displayed the resulting FROG trace 
and resized and conditioned the FROG trace for the phase retrieval process. 
The computer also was the user interface, displaying the retrieved pulse and 
the raw video image. VideoFROG used the PCGP phase retrieval process, 
which generally converges well. Consequently, no provision was made to 
detect stagnation or other faults in the phase retrieval process. 



SUMMARY OF THE PREFERRED EMBODIMENTS 
[0029] A preferred implementation may include a real-time pulse 

measurement device that has direct, integrated feedback that measures how 
well the device is retrieving pulses and tracking changes in the pulse train. 
This feedback may be in the form of the FROG trace error, the display of the 
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measured and retrieved FROG trace, or another display or algorithm that 
can be used to provide information about the quality of the measurement. 
[0030] According to one aspect, the present invention may provide a 
method for real-time measurement of ultrashort laser pulses comprising 
recording in a computer measured frequency resolved optical gating (FROG) 
trace data, the FROG trace data generated by processing a pulse in a FROG 
device. The method continues by processing the measured FROG trace to 
perform real time phase retrieval and to generate in real time a retrieved 
pulse from the measured FROG trace. The retrieved pulse is displayed. A 
feedback parameter is generated in real time to provide information on the 
functioning of the real time phase retrieval. 

[0031] According to another aspect, the present invention may provide 

a method of performing real time phase retrieval processing of frequency 
resolved optical gating (FROG) traces. The method comprises receiving as 
input a measured FROG trace data set, the FROG trace data set generated 
by processing a pulse in a FROG device and processing the measured FROG 
trace data set to perform real time phase retrieval and generating in real 
time a retrieved pulse from the measured FROG trace. Displays of the 
retrieved pulse are preferably generated at a rate of three Hertz or faster. A 
feedback parameter provides information characterizing the real time phase 
retrieval. 
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[0032] A method of performing real time phase retrieval processing of 
frequency resolved optical gating (FROG) traces, the method comprising 
receiving as input a measured FROG trace data set created by processing a 
pulse in a FROG device. The method includes processing the measured 
FROG trace data set to perform real time phase retrieval and generating in 
real time a retrieved pulse from the measured FROG trace. Real time 
displays are generated of the retrieved pulse. The method also includes 
generating in real time a feedback parameter providing information 
characterizing the real time phase retrieval and performing a control 
operation in response to the feedback parameter. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0033] The accompanying drawings, which are incorporated into and 

form a part of the specification, illustrate several embodiments of the present 
invention and, together with the description, serve to explain the principles of 
the invention. The drawings are only for the purpose of illustrating a 
preferred embodiment of the invention and are not to be construed as 
limiting the invention. 

[0034] FIG. 1 illustrates that measuring the spectrogram of a pulse is 

easier than measuring its sonogram. A spectrogram is a spectrally resolved 
autocorrelation. In this figure, the optical Kerr-effect is used (polarization- 
gate) as the non-linearity. 
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[0035] FIG. 2 shows schematically an SHG FROG device. The SHG 

signal from the autocorrelation is spectrally resolved. SHG FROG is very 
simple and sensitive, but has a direction-of-time ambiguity. For example, if 
the pulse has chirp, only the magnitude of the chirp is determined — the sign 
of the chirp remains unknown. 

[0036] FIG. 3 illustrates phase retrieval processes for the inversion of 

FROG spectrograms that start with an initial guess for the pulse to generate 
an initial E S i g (t,T). A I'D Fourier transform generates the FROG trace. The 
process continues by replacing the magnitude of the calculated FROG trace 
with the square root of the measured FROG trace. The process takes the 
inverse Fourier transform with respect to © to produce the new signal field 
and generate a new start value for E(t). Interestingly, it is only the process 
that produces E(t) from E S i g (t,x) that differentiates all the FROG 
implementations. 

[0037] FIG. 4 is a schematic illustration of a FROG system in 

accordance with aspects of the present invention. 

[0038] FIG. 5 illustrates the display and preprocessing of raw FROG 

data, along with the iterative phase retrieval process used to characterize a 
high-speed laser pulse. According to preferred aspects of the present 
invention, the illustrated FROG system provides an indication of feedback 
that can be used to identify an error condition and one of the responses to an 
error condition, restarting phase retrieval using a new form of initial guess. 
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[0039] FIG. 6 illustrates the optical path for an exemplary polarization 

gate FROG device that could be used in the FROG system of FIG. 4. 
[0040] FIG. 7 illustrates the optical path for an exemplary second 

harmonic generation (SHG) FROG device that could be used in the FROG 
system of FIG. 4. 



DESCRIPTION OF THE PREFERRED EMBODIMENTS 

[0041] Preferred embodiments of the present invention provide a real- 

time ultrafast laser pulse measurement method and apparatus, which 
overcomes the low speed of pulse measurement using frequency-resolved 
optical gating (FROG). Particularly preferred implementations identify 
stagnation, loss of tracking or other errors in a real-time pulse measurement 
system based on frequency-resolved optical gating. Some embodiments of the 
present invention preferably condition the spectrogram or FROG trace for 
better and more accurate phase retrievals. Other aspects of the present 
invention preferably calculate parameters that indicate the quality of the 
FROG measurement. 

[0042] The PCGP phase retrieval technique discussed in the 

background generally converges well, but it is possible for the PCGP phase 
retrieval technique to stagnate, lose tracking or for other errors to develop. 
These errors may result from noise, unremovable background or poor 
calibration and may be particularly problematic for inexperienced users of 
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the FROG system. It is also possible to use other phase retrieval strategies 
to provide real-time FROG, which are more prone to stagnation and errors 
than PCGP phase retrieval. Consequently it is preferable to provide a real 
time pulse measurement system that provides feedback to a user to allow 
identification of stagnation or other types of errors. 

[0043] Preferred implementations of this real-time pulse measurement 

system can measure ultrafast laser pulses in real-time at rates on the order 
of a few Hertz, such as 3 Hz, or at even faster rates such as 20 Hz or faster. 
The method and apparatus preferably integrates a frequency-resolved optical 
gating (FROG) device with a desired phase retrieval process in a seamless 
way. The data signal, or FROG trace, is collected using a suitable data 
acquisition device such as a video camera. The video signal is digitized using 
a frame grabber and is displayed on the computer screen. This signal, the 
raw FROG trace, provides the user with an indication as to the performance 
of the FROG device itself. Because the raw FROG trace is a spectrogram of 
the pulse - a display of intensity vs. frequency and time, it also provides some 
indication as to the duration and frequency content of the pulse. To obtain 
the actual intensity and phase of the pulse, a two-dimensional phase retrieval 
algorithm is applied to the FROG trace. 

[0044] To use the phase retrieval algorithm the raw FROG trace is 

resized to establish the correct relationship between the time and the 
frequency axis such that the frequency span is equal to the inverse of the 
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time spacing and the time span is equal to the inverse of the frequency 
spacing, which is the general case for Fourier transforms. 
[0045] After the FROG trace is resized, it is subjected to the FROG 
processing. To improve the speed of the tracking of the FROG processing, the 
pulse retrieved from the previous FROG trace is used as the initial guess for 
the new FROG trace. In this way, the FROG algorithm can be thought of as 
an adaptive digital filter. While the PCGP algorithm has usually been used 
in real-time FROG measurements, other FROG retrieval algorithms can be 
used as well. 

[0046] Because large changes can occur in the pulse even at the 

retrieval rates possible using this method and apparatus, the phase retrieval 
process can stagnate and lose tracking of the pulse. Indeed, this can even 
happen if the pulse does not change. Previous real-time pulse measurement 
systems provide no feedback to the user to indicate a problem. In order to 
provide feedback to both the program and the user, preferred embodiments of 
the present invention provide within the preferred program implementation a 
calculated and/or displayed metric to indicate how well the iterative FROG 
process is tracking the pulse. One preferred choice and typically the simplest 
metric to use is the FROG trace error. High errors indicate a loss of tracking; 
that is, the retrieved pulse is no longer close to the true pulse. Consequently, 
the FROG phase retrieval process is preferably reset. A particularly 
preferred strategy for resetting the FROG phase retrieval process is to adopt 



19 



a Gaussian pulse that has random phase as the initial guess in the process 
resumed after the loss of tracking or stagnation. 
[0047] As previously discussed, the FROG trace error metric is 
proportional to the distance measurement between the FROG trace 
constructed from the retrieved pulse (retrieved FROG trace) and the 
measured FROG trace that has been resized. Typically, if the measured 
FROG trace is close to the retrieved FROG trace, then the distance is small, 
and the phase retrieval process is tracking properly. If the distance becomes 
large, then the phase retrieval process may be stagnated or may no longer be 
tracking. The distance measurement is usually the Euclidean norm of the 
difference between the measured and the retrieved FROG trace, divided by 
N, where N is the square root of the number elements in the FROG trace. 
Other distance related measurements can be used as the FROG trace error 
metric such as the sum of the absolute value of the differences (the h norm of 
the difference) or the per element average of the absolute value of the 
differences. 

[0048] While small distances (low FROG trace errors) between the 

measured and retrieved FROG traces usually indicate good retrievals, they 
are not always a complete indicator of stagnation or a loss of tracking. These 
quantities are subject to interpretation. For example, retrievals from large 
FROG traces always have smaller errors than comparable retrievals from 
smaller FROG traces. This usually occurs because larger FROG traces have 



20 



more elements with small values than smaller FROG traces. Also, if the 
FROG trace is very complex and has large amounts of structure, then the 
FROG trace error may be higher than for a smooth, simple FROG trace even 
when both FROG traces are nearly the same size. For example, an 
acceptable FROG trace error is in the range of 0.5% - 1.5% for second 
harmonic generation (SHG) FROG retrievals on a 64 x 64 grid. But errors as 
high as 3.5% can still yield acceptable results. If the grid size is increased to 
128 x 128 then the acceptable error drops by a factor of 2. For polarization 
gate FROG, acceptable FROG trace errors are on the order of 2%. Again, for 
larger grid sizes, this value drops. Thus, the FROG trace error alone 
provides a guide to the quality of the retrieval, but is subject to 
interpretation. 

[0049] Another way to determine if the real-time FROG device is 

measuring pulses correctly is to visually compare the measured FROG trace 
with the retrieved FROG trace. When the two FROG traces are compared 
side-by-side, they should appear very similar. Thus, if the measured and 
retrieved FROG traces appear similar and the FROG trace error is low, then 
the retrieval is likely of good quality. If they have the same general shape 
and the FROG trace error is a little high, then the retrieval is probably good. 
If the two traces appear very different, then the retrieval is almost certainly 
poor. 
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[0050] While it might be possible to automate such feedback and 

analysis, automation is quite difficult in practice. A carefully constructed 
neural net might be able to make a correct decision about the quality of the 
retrieval, but having the displays easily available for the user is presently 
preferred as a means of providing an indication of the quality of the retrieval 
and, therefore, the quality of the measured pulse. In addition, the displays 
provide the user with enough information to tell if the algorithm is no longer 
tracking the pulse. 

[0051] A preferred implementation of a FROG system is shown in 
FIG. 4. A short pulse 10, such as an ultrashort laser pulse, is sent into a 
frequency resolved optical gating (FROG) device 12. The spectrogram output 
of the FROG device is recorded using a video camera 14. The video output 
from the video camera 14 is digitized by a frame grabber 16, and sent into the 
computer 18, where the data are analyzed and displayed for the user. FROG 
processing and phase retrieval are preferably accomplished in computer 18, 
which preferably implements the PCGP technique described in U.S. patent 
No. 6,219,142, which is incorporated by reference in its entirety for all of its 
teachings related to FROG phase retrieval and FROG apparatus. 
[0052] It should be appreciated that the various preprocessing and 

feedback strategies described here are preferably implemented in conjunction 
with the PCGP phase retrieval technique in a FROG system. It should 
further be appreciated that the preprocessing and feedback strategies 
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described here can also be implemented advantageously in FROG 
implementations that do not use the PCGP technique. In fact, aspects of the 
present invention are expected to have even greater benefits in FROG 
implementations that do not use PCGP techniques. 
[0053] The analysis process is shown in FIG. 5. FIG. 5 shows the 

display and preprocessing of raw FROG data, along with the iterative phase 
retrieval process used to characterize the pulse. According to preferred 
aspects of the present invention, the illustrated FROG system provides an 
indication of feedback that can be used to identify an error condition and one 
of the preferred responses to an error condition, restarting phase retrieval 
using a new form of initial guess. Data obtained by the frame grabber 16 
(also shown in FIG. 4) are provided to the computer 18 (FIG. 4) for processing 
and display. The computer both displays 22 (FIG. 5) and conditions 24 the 
raw data for the iterative two-dimensional phase retrieval process 26. Prior 
to sending the new conditioned data 24 to the iterative retrieval process 26, 
the results from the previous iteration are extracted and displayed 28. 
[0054] A feedback parameter is calculated 32 and fed back 34 to control 

the phase retrieval process 26. One possible option for control is to reset the 
phase retrieval process by using a new initial guess 36, for example using a 
Gaussian pulse that has random phase as the initial guess 36 in the process 
resumed after the loss of tracking or stagnation. Another strategy is for the 
program to alter the phase retrieval process, for example by introducing a 
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small, random or predetermined, perturbation into the then-current 
parameters of the phase retrieval process and starting another iteration 
using the perturbed information. 

[0055] Another sort of action that might be taken in response to the 

feedback signal is for the system to generate an alarm or alert error 
indication. Generally a visual alarm or alert is preferred, although an 
audible alarm or alert may alternately be used or may be used in addition to 
a visual indication. 

[0056] The feedback parameter 34 for control of the algorithm 26 can 

be one or more of several different parameters. It could be the FROG trace 
error where if this error is too great, the phase retrieval process has lost track 
and must be reset. As discussed above, a FROG trace error metric used as 
the feedback parameter may be proportional to the distance measured 
between the FROG trace constructed from the retrieved pulse and the 
measured FROG trace that has been resized. Any distance related 
measurement could be used. Usually the Euclidean norm, divided by the 
square root of the number of elements in the FROG trace, is used. 
[0057] A different feedback parameter 34 could be a comparison of the 

conditioned measured FROG trace with the trace constructed from the 
retrieved pulse, which is called the retrieved FROG trace. As discussed 
above, this comparison is preferably provided directly to the display of the 
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computer, where a skilled user can efficiently determine the status of the 
phase retrieval process. 

[0058] The feedback parameter 34 could also be a metric measuring the 

amount of background removed by background subtraction. For this 
background-related feedback parameter, the parameter preferably identifies 
conditions where too much background has been removed and conditions 
where the background in the signal is too high. This background-related 
feedback parameter is discussed in greater detail below. 
[0059] FIG. 6 illustrates the optical path of a polarization gate FROG 

device 12 (FIG. 4) that might be used in the FROG system of FIG. 4. An 
ultrafast laser pulse 40 to be characterized is provided to the FROG device 
and is split by a beam splitter 42. A portion of the input pulse, 40 passes 
through the beam splitter 42 and is directed by mirrors including 
retroreflector 44 through a first polarizer 46, which polarizes the probe 
portion 48 of the pulse. A spherical lens 50 focuses the polarized portion 48 of 
the pulse onto a non-linear medium 52. Typically the non-linear medium 52 
is quartz. Mirror 54 directs the portion of the input pulse 40 that is split off 
and reflected from the beam splitter along an optical path so that lens 50 
focuses this gate portion of the pulse 56 onto the non-linear medium 52. 
[0060] Two resulting pulses 48, 56 form the probe and time delayed 

gate pulses that are combined by the non-linear medium 52. The gate pulse 
56 produces a birefringence in the nonlinear material 52, which causes the 
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polarization of the probe pulse 48 to rotate slightly allowing the signal pulse 
62 to pass through polarizer 60 that is crossed relative to the input polarizer 
46. A spatial filter 58 and the polarizer 60 remove unwanted portions of the 
initial pulses to produce the signal pulse 62 shown conceptually in the inset 
64. The inset 64 shows the temporal relation between the probe, signal and 
gate pulses, as well. Lens 66 and mirror 68 direct the signal pulse 62 into the 
spectrometer and detector 70 that produces the two-dimensional spectrogram 
of the signal pulse. Video camera 14 (FIG. 4) and frame grabber 16 (FIG. 4) 
capture the spectrogram for phase retrieval processing as discussed above. 
[0061] FIG. 7 illustrates the optical path of a second harmonic 

generation (SHG) FROG device 12 (FIG. 4) that might be used in the FROG 
system of FIG. 4. An input pulse 80 enters the second harmonic generation 
FROG device and is split into two replicas by a beam splitter 82. One pulse, 
which is deemed the gate pulse 84, is retro reflected by two mirrors 86. The 
other replica, the probe pulse 88 is sent to a movable delay stage 90. A 
mirror 92 directs the probe pulse 88 through a spherical lens 94 that focuses 
both the probe pulse 88 and the gate pulse 84 into a nonlinear crystal 96 
suitable for second harmonic generation. The probe pulse 88 and the gate 
pulse 84 sum in the nonlinear crystal 96 to produce the signal pulse 98. A 
spherical lens 100 images the signal pulse 98 onto the slits of a spectrometer 
102. An optional mirror 104 helps to direct the signal pulse 98 into the 
spectrometer 102. An iris 106 prevents any stray second harmonic light 
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generated by the probe and gate pulses from entering the spectrometer 102 
and corrupting the signal pulse 98. The detector records the spectra as a 
function of delay to produce the FROG trace, which is subjected to phase 
retrieval processing as discussed above. 

[0062] The preceding discussion has emphasized errors associated with 

stagnation or loss of tracking and the use of feedback to address these errors. 
FROG trace retrievals are also subject to other issues that affect the retrieval 
quality. An important problem with FROG retrievals is the deleterious effect 
noise in the wings of the measured FROG trace can have on the retrieved 
pulse. Such noise can cause artifacts in the retrieved pulse and can lower the 
dynamic range of the measurement. While these artifacts are most 
noticeable in the wings of the pulse, where the pulse intensity is small, these 
background-induced artifacts occur randomly anywhere on the pulse. 
[0063] Several strategies can be used to mitigate background effects. 
One way is to tailor the camera response such that the camera is very 
responsive to low light levels and less responsive to high light levels. Such a 
response is already available on video cameras in the form of the gamma 
compensation or correction function. Video cameras provide gamma 
correction to compensate for the non-linear response of phosphors in cathode- 
ray tubes (CRTs). Gamma correction is an analog processing that is 
implemented in certain aspects of the invention to condition the FROG trace 
for processing in a manner that improves the dynamic range of the data. 



27 



[0064] When the gamma correction is engaged, the recorded intensity is 

approximately the light intensity raised to a fractional power, R = F, where R 
is the response, I is the light intensity on the camera, and y, or "gamma" is 
the gamma correction. The true light intensity on the camera is then I=R 1/y . 
When gamma is less than 1, the gamma correction has the effect of 
increasing the range of the signal and suppressing any low level noise, which 
usually appears in the wings of the FROG trace. The industry standard 
gamma is 0.45, close to 0.5. Indeed, early real-time pulse retrieval devices 
used this to provide an analog calculation of the square root of the intensity, 
which is required by the phase retrieval algorithm. 

[0065] Actual gamma corrections on cameras are not 0.45, however. 

They are usually closer to 0.6 - 0.7 even when they are specified to be 0.45. 
Consequently, one cannot assume the gamma is exactly 0.45 or 0.5. The true 
gamma is preferably measured and the exact inverse of the gamma is most 
preferably applied to the recorded signal. However, even though the gamma 
correction cannot be used as an estimate of the square root of the intensity, it 
can be used to suppress the effect of noise in the wings of the FROG trace. 
Thus, implementations of the present invention that utilize this aspect of the 
invention preferably utilize a gamma correction in the video camera (14 in 
FIG. 4). The level of gamma correction is measured and the actual gamma 
correction is preferably stored for future uses by the preferred program 
implementations of this aspect of the present invention. The gamma 
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correction preferably is reversed in a calculation within the computer to 
produce the original intensity distribution of the FROG trace before phase 
retrieval processing. Preferably, the use of gamma correction is selectable by 
the user through the user interface of the computer. 

[0066] Another way the effect of noise can be mitigated is to filter the 
FROG trace. Generally this filtering is done in the computer 18 prior to 
phase retrieval processing, but it is possible to implement the filtering in 
different ways including within the frame grabber. One useful type of filter 
that can be used to remove speckle noise in the wings of the FROG trace is 
the median filter. A median filter works by replacing the center of a small 
sub-region of the image by the median of the values in the sub-region. This is 
accomplished by ordering the values in a small (N x N) sub -region of an 
image or matrix from lowest to highest (or vice versa). If N is an odd number, 
the total number of points is odd, and the middle value in the ordered 
sequence is by definition the median value. The sub-region is moved over the 
entire image so that every point of the image is filtered. If the center point of 
the sub-region is on the edge of the matrix or image, then wrapping the 
image - using points from opposite sides of the matrix, can fill in points in 
the sub-region. Alternatively, fixed values, such as zero, can be added to the 
sides of the image to fill the sub-region when it is on the edge of the image. 
[0067] Another type of filter that can be used is a rank filter. Rather 
than choose the median, or middle value, of an N x N sub-region of the image, 
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a value with lower or higher rank is chosen, such as the lowest value, the 
next to lowest value, the highest value, etc. In addition, a modification of the 
rank filter determines if the center point of the sub-region is closer to the 
lowest or highest value in the sub-region. If the point is closer to the highest 
value, it is replaced by the higher value. Alternatively, if it is closer to the 
lowest value, it is replaced by the lowest value. 

[0068] For purposes of filtering FROG traces, it is often best to make 

the rank filter adaptive by applying the rank filter for only certain values of 
the FROG trace. Most often, the FROG trace needs to be filtered by a rank 
filter when the values are low. When the rank filter is made adaptive in this 
manner, it requires two numbers to be specified: the rank and the threshold. 
If the sum of the values in the sub-region is greater than the threshold, 
nothing is done. Conversely, if the sum of the values in the sub-region is less 
than the threshold, the center value of the sub-region is replaced by the value 
with the appropriate rank. Preferably, the use of any of these filter functions 
is selectable by the user through the user interface of the computer. 
[0069] This filter is not only useful for FROG traces, but it is also 

useful for situations where noise in the wings of an image interferes with a 
calculation. Most notably, noise in the wings of an image of a laser beam 
interferes with the calculation of the second moment of the laser beam, which 
is often used to specify the diameter of the laser beam. Any noise or speckle 
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in the wings of the beam image (profile) may cause the second moment 
calculation to have an error of several percent. 

[0070] Because background has such deleterious effects on the retrieval 

of the pulse from its FROG trace, users tend to remove too much background, 
causing errors in the retrieved pulse. Thus, another important tool in 
improving the quality of FROG trace retrieval is the ability to determine if 
too much background has been removed from the measured FROG trace. 
[0071] Typically, the FROG trace does not go directly to zero around its 

perimeter; it asymptotically approaches zero. If the data acquisition were 
perfect, then these wings would be visible and noise free. However, data 
acquisition devices have detection limits and the wings of the FROG trace 
may fall under these limits. Background removal may further exacerbate 
this problem to the point where significant portions of the FROG trace are 
removed and, as a result, the retrieval can appear to be non-physical. 
[0072] When significant amounts of background are missing from a 

FROG trace, the remaining FROG trace can be thought of as a cap on the 
true FROG trace. That is, if the measured FROG trace is superimposed on 
the true FROG trace, then while the peak(s) of the measured and true FROG 
trace would coincide, the measured FROG trace would not reach the wings of 
the actual FROG trace and the edge of the measured FROG trace would sit 
above the bottom of the true FROG trace by some amount p. This can be 
visualized by thinking of a mountain surrounded by an ocean where ocean 
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level represents the detection limit of the data acquisition device. The 
portion of the mountain above the ocean represents the measured FROG 
trace while the entire mountain represents the true FROG trace. The 
distance between the surface of the ocean and the base of the mountain is p. 
[0073] Typically, when the phase retrieval process is applied to the 

measured FROG trace, the retrieved pulse is such that the retrieved FROG 
trace more closely approximates the true FROG trace than the measured 
FROG trace. That is, the retrieval process fills in some of the wings clipped 
from the measured FROG trace. If the measured FROG trace is then 
normalized to the retrieved FROG trace such that the two heights are the 
same, the FROG trace error will be abnormally high because portions of the 
measured FROG trace are zero when they should actually be non-zero. 
[0074] However, if p is added (e.g., in computer 18 of FIG. 4) to the 

measured FROG trace before the height of the measured FROG trace is 
normalized to the height of the retrieved FROG trace, and the error is only 
calculated using points where the measured FROG trace is non-zero, the 
FROG trace error will drop. Indeed, an estimate for p can be found by finding 
the offset that produces the lowest FROG trace error. 
[0075] Feedback to the user is not just useful as an indication of 

measurement and system function. Feedback can also be used to adjust the 
ultrafast laser source for optimum performance. To facilitate the feedback, 
information from the retrieved pulse may have to be conditioned. For 
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example, in order to amplify ultrafast laser pulses, the ultrafast laser pulses 
is preferably lengthened in time, which is accomplished by chirping the pulse 
using a pulse stretcher. After amplification, a pulse compressor recompresses 
the pulses. Because the stretcher and compressor must be closely matched, 
they are usually similar in design. Typically, the stretcher and compressor 
are made by dispersing a beam off a grating using a focusing element to 
collimate the light to another focusing element that focuses the beam back 
onto another grating before being output. 

[0076] The distance between the gratings determines the group velocity 

dispersion, or linear chirp. Third order dispersion is a function of the 
difference between the gratings angle used in the system. Adding prisms to 
the system can compensate for fourth order dispersion. From the frequency 
domain phase, the amount of group velocity dispersion, as well as all higher 
order dispersion terms can be found and displayed for the user or the 
information can be used to automatically adjust the stretcher-compressor 
system. The information could also be sent to the factory, electronically or 
otherwise, to assist with maintenance or as a diagnostic of the ultrafast laser 
amplifier. 

[0077] Real time as it is used in this application means a rate of at 

least one Hertz. Rates of 3 Hz or faster are more preferred real time 
implementations to provide a useful rate of update and rates of 10 Hz or 
better are further preferred. As discussed above, rates of 20 Hz or faster are 
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achievable in particular preferred embodiments of the present invention and 
that update rate is a still further preferred update rate at the present time. 
[0078] Useful information is not just limited to the frequency domain, 

or spectral phase of the pulse. Time domain phase can be used to determine 
spectral loss of the pulse due to clipping of the beam within the stretcher- 
compressor. 

[0079] This device may be used for measuring light pulses in 

telecommunications systems as well. 

[0080] Although the present invention has been described in detail 
with reference only to the presently preferred embodiments, those of ordinary 
skill in the art will appreciate that various modifications can be made 
without departing from the invention. Accordingly, the invention is not to be 
limited to any of the described embodiments thereof but is instead defined by 
the following claims. 



34 



